Negative circuits and sustained oscillations 
in asynchronous automata networks 



Adrien Richard 



Laboratoire I3S, UMR 6070 CNRS & Universite de Nice-Sophia Antipolis, 
2000 route des Lucioles, 06903 Sophia Antipolis, France. 



e-mail: richard@i3s.unice.fr 

telephone: +33 4 92 94 27 51 
fax: +33 4 92 94 28 98 



1 



Abstract: The biologist Rene Thomas conjectured, twenty years ago, that the presence 
of a negative feedback circuit in the interaction graph of a dynamical system is a necessary 
condition for this system to produce sustained oscillations. In this paper, we state and 
prove this conjecture for asynchronous automata networks, a class of discrete dynamical 
systems extensively used to model the behaviors of gene networks. As a corollary, we obtain 
the following fixed point theorem: given a product X of n finite intervals of integers, and a 
map F from X to itself, if the interaction graph associated with F has no negative circuit, 
then F has at least one fixed point. 
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1 Introduction 

We are interested in a class of discrete dynamical systems used to model gene networks. 
The biological context is the following. Gene networks are often described by Biologists 
under the form of interaction graphs. These are directed graphs where vertices correspond 
to genes and where arcs are labelled with a sign: a positive (negative) arc from a gene j to 
a gene i means that the protein encoded by the gene j activates (represses) the synthesis 
of the protein encoded by the gene i. These very coarse descriptions of gene networks are 
then taken as a basis to design much more complex dynamical models that describe the 
temporal evolution of the concentration of the encoded proteins [TJ. Unfortunately, these 
models require, in most cases, unavailable informations on the strength of the interactions. 
In this context, a difficult and interesting question is: which dynamical properties of a gene 
network can be deduced from its interaction graph ? 

The biologist Rene Thomas stated two well known conjectures that partially answer 
this question. These conjectures can be informally stated as follows \25\ [9]: 

1. The presence of a positive circuit in the interaction graph of a network (i.e. a circuit 
with an even number of negative arcs) is a necessary condition for the presence of 
multiple stable states in the dynamics of the network. 

2. The presence of a negative circuit in the interaction graph of a network (i.e. a circuit 
with an odd number of negative arcs) is a necessary condition for the presence of 
sustained oscillations in the dynamics of the network. 

It is worth noting that multistationarity and sustained oscillations are, from a biological 
point of view, important dynamical properties often related to differentiation processes and 
homeostasis phenomena respectively [25 | [26 ] [25]. 
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The first conjecture has been formally stated and proved by several authors in contin- 
uous frameworks |101 [6j [2T| [H [22], [23] , in which the concentration of each protein evolves 
continuously, generally following an ordinary differential equation system. The first con- 
jecture has been more recently stated and proved in discrete frameworks pi [2] [TBI HT] 
in which the concentration level of each protein evolves inside a finite interval of integers, 
which is {0, 1} in the Boolean case. Studies of the second conjecture are fewer: a Boolean 
version of the second conjecture has been stated and proved by Remy, Ruet and Thieffry 
[TT] . and there are only partial results in the continuous case [6, 2 lj . 

In this paper, we state and prove second Thomas' conjecture for asynchronous automata 
networks (Theorem 1). Our interest for these discrete dynamical systems comes from the 
fact that they have been proposed by Thomas as model for the dynamics of gene networks 
more than thirty years ago [2U [261 113 EE] • They are still extensively used because of the 
qualitative nature of most reliable experimental data, and the fact that the sigmoidal shape 
of genetic regulations leads to a natural discretization of concentrations [5j [191 EH EQl [H] . 

The discrete version of Thomas' conjecture we establish generalizes in several ways the 
one established by Remy, Ruet and Thieffry [TTJ in the Boolean case: both the discrete 
dynamical framework and the considered class of sustained oscillations are more general. 
Furthermore, the class of sustained oscillations we consider allows us to obtain, as an 
immediate consequence, the fixed point theorem mentioned in the abstract (Corollary 1). 

The paper is organized as follows. Section 2 presents definitions related to asynchronous 
automata networks. In Section 3, second Thomas' conjecture is stated and proved for these 
networks. In Section 4, we establish a variant of second Thomas' conjecture more suited to 
the modeling of gene networks. Counter examples to natural extension of the established 
results are given in Section 5. 

4 



2 Definitions 

We consider a network of n interacting automata, denoted from 1 to n. The set of possible 
states for automaton i is a finite intervals of integers X{ of cardinality at least two. The set 
of possible states for the network is the Cartesian product X = JJiLi ^i- The dynamics of 
the network is then described according to a map F : X — > X, 

x = (xi, . . . ,x n ) G X i > F(x) = (/i(x), . . . ,f n (x)) G X, 

with which we associate the maps F{ : X — > X defined by 

Fi(x) = (xi, . . . ,Xi_i,/i(x),x i+ i, . . . ,x n ) (i = l,...,n). 

More precisely, given an initial point x° G X and a map ip from N to {1, . . . , n}, the dynam- 
ics of the network is described by the following recurrence, that we call the asynchronous 
iteration of F induced by the strategy ip from initial point x°: 

x t+1 = F^x*) (t = 0,1,2...)- (1) 

Generally, one only considers the asynchronous iterations induced by pseudo-periodic strate- 
gies, i.e. strategies tp such that |</? _1 (z)| = oo for i = 1, . . . , n [171 [3]. 

In this paper, we will study the asynchronous iterations of F through a directed graph 
on X called the asynchronous state transition graph of F. Before defining this graph, let 
us set, for all x G X, 

I F (x) = {ie{l,...,n} \ fi(x) / x^. 
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Definition 1 The asynchronous state transition graph of F, denoted T(F), is the directed 
graph whose set of vertices is X and whose set of arcs is 

{(x,Fi(x))\x £l, ieI F (x)}. 

Remark 1 is the number of successors of x in T(F), and |/p(a;)| = if and only if 

x is a fixed point of F. Also, T(F) has no arc from a vertex to itself, and in the following, 
we assume, by convention, that T(F) has a path of length zero from each vertex to itself. 

The relation between T(F) and the asynchronous iterations of F is clear: there is a path 
from x to y in T(F) if and only if there exists a strategy ip such that the asynchronous 
iteration of F induced by ip from x reaches y. 

In this context, the fixed points of F are of particular interest: they correspond to the 
stable states of the system. More precisely, if <p is a pseudo-periodic strategy, then the 
asynchronous iteration (pTJ) stabilizes on a point £ (i.e. there exists t such that x l = x t+l = 
£) if and only if £ is a fixed point of F. In the following definition, we introduce a notion 
of an attractor, which extends in a natural way the one of a stable state. 

Definition 2 A trap domain of T(F) is a non-empty subset D C X such that for every 
arc (x,y) ofT(F), if x G D then y € D. An attractor ofT(F) is a smallest trap domain 
with respect to the inclusion. A cyclic attractor is an attractor of cardinality at least two. 

Remark 2 One has the three following basic properties: (1) x is a fixed point of F if and 
only if {x} is an attractor of T(F); (2) attractors perform an attraction in the weak sense 
that, from any state, there always exists a path leading to one of them; (3) if x and y 
belong to the same attractor, then there exists a path from x to y. 
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The third point highlights the fact that inside a cyclic attractor, each state has at least 
one successor. So, when the network is inside a cyclic attractor, it cannot reach a fixed 
point, and thus, it describes sustained oscillations. More precisely, if x° belongs to a cyclic 
attractor A, then for all pseudo-periodic strategy if, the asynchronous iteration of F induced 
by (p from x° never leaves A and never stabilizes, and since A is finite, it necessarily 
describes sustained oscillations. In the following, we are interested in the relationships 
between sustained oscillations produced by cyclic attractors and the negative circuits of 
the interaction graph of the network. 

An interaction graph is here defined to be a directed graph whose set of vertices is 
{1, . . . ,n} and where each arc is provided with a sign. Formally, each arc is characterized 
by a triple (j, s, i) where j (i) is the initial (final) vertex, and where s G {—1, 1} is the sign 
of the arc. An interaction graph can then have both a positive and a negative arc from 
one vertex to another. 

In the following definition, we attach to F an interaction graph G(F) that is nothing 
but the interaction graph of the network whose dynamics is described by the asynchronous 
iterations of F. 

Definition 3 The interaction graph of F, denoted G(F), is the interaction graph that 
contains a positive ( negative ) arc from j to i if there exists x£l with Xj + 1 £ Xj such that 

fi{%li ■ ■ ■ i Xj + 1, . . . , X n ) fi{x\, . . . , Xj , . . . , X n ) 

is positive (negative). 

Remark 3 G(F) has at least one arc from j to i if and only if /j depends on Xj. 
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Definition 4 A path of G(F) of length r > 1 is a sequence of r arcs of G(F), say 
(ji, si, h), (j2, 82, 12), ■ ■ ■ , (jr, s r , i r ), such that i g = j q+ \ for all 1 < q < r. Such a path 
is a path from j\ to i r of sign s = YY q =i s q- It ^ a circuit if i r = j\ and it is an elemen- 
tary circuit if, in addition, the vertices i q are mutually distinct. 

Remark 4 If G[F) has a negative circuit, then it has an elementary negative circuit (this 
is false for positive circuits). So, in order to prove that G(F) has an elementary negative 
circuit, it is sufficient to prove that G(F) has a negative circuit. 

Example 1 n = 2, X = {0, 1, 2} 2 and F is defined by the following table: 



X 


(0,0) 


(0,1) 


(0,2) 


(1,0) 


(1,1) 


(1,2) 


(2,0) 


(2,1) 


(2,2) 


F(x) 


(2,0) 


(1,0) 


(0,2) 


(2,0) 


(0,0) 


(0,1) 


(2,1) 


(0,1) 


(0,1) 



The asynchronous state transition graph and the interaction graph of F are as follows: 
F(F) 



(0,2) < (1,2) (2,2) 



G(F) 



(0,1) ^±(1,1) (2,1) 

A 



(0,0) (1,0) >(2,0) 




-1 



We see that T(F) has two attractors: the stable state (0, 2) and the cyclic attractor {0, 1, 2}x 
{0,1}. We also see that G(F) has two elementary positive circuits, and two elementary 
negative circuits. 



3 Main result 

In this section, we prove the following discrete version of second Thomas' conjecture: 

Theorem 1 IfT(F) has a cyclic attractor, then G(F) has a negative circuit. 

Remark 5 This theorem has been proved by Remy, Ruet and Thieffry [11] in the Boolean 
case (i.e. when X is the n-dimensional hypercube {0, l} n ) and under the rather strong 
hypothesis that T(F) contains a stable cycle, that is, a cyclic attractor A in which each 
state has a unique successor (i.e. = 1 for all x E A). 

Before proving Theorem 1, let us point out that it has, as immediate consequence, the 
following fixed point theorem (which can not be deduced, in the Boolean case, from the 
theorem of Remy, Ruet and Thieffry mentioned above): 

Corollary 1 If G(F) has a no negative circuit, then F has at least one fixed point. 

Proof — Indeed, if F has no fixed point, then T(F) has clearly at least one cyclic attractor, 
and following Theorem 1, G(F) has a negative circuit. □ 

Remark 6 In [17, Chapter 13] (see also [3]), Robert prove the following convergence result: 
if G(F) has no circuit, then F has a unique fixed point £, and, for all initial point x° and 
for all pseudo-periodic strategy <p, the asynchronous iteration of F induced by <p from x° 
reaches the fixed point £. Prom Theorem 1 and the second point of Remark 2, one obtains 
a convergence result that has a weaker conclusion under a weaker condition: if G(F) has 
no negative circuit, then F has at least one fixed point, and for all initial point x°, there 
exists a strategy <p for which the asynchronous iteration (pQ) reaches a fixed point of F. 
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The proof of Theorem 1 needs few additional definitions and notations. Let Q and Q' 
be two interaction graphs with arc-set E and E' respectively. We say that Q is a subgraph 
of Q' if E C E' . We denote by Q U Q' the interaction graph whose set of arcs is E U E' . 
Next, for all x € X, we set 

/■(x) = sign(/i(x) - Xj) (i = 1, . . . , n), 

where sign is the usual sign function (sign(a) = a/\a\ for all and sign(O) = 0). The 

main tool used in the proof of Theorem 1 is the following notion of local interaction graph: 

Definition 5 For all x € X, we denote by ^p(x) the interaction graph that contains an 
arc from j to i of sign s G { — 1, 1} if 

fi(.x) + fl(Fj(x)) and s = /j(x)/^(x)) 
Lemma 1 For all x € X, 5fp(x) is a subgraph of G(F). 

Proof — Let x G X, and suppose that £^f(x) has an arc from j to i of sign s. For all 
integer p, we set 

X P — {x\ , ••• , Xj — 1 , Xj + p, Xj + \ , . . . , x n ) . 

Case j 7^ i. By definition, /j(x) ^ 0. We suppose /j(x) > the other case being 
similar. Setting q = fj{x) — Xj, we have q > and x q = Fj(x). So s = fl(x q ) and 
/Kz) = /^O) ^ Consider the smallest < p < 9 such that fl(x p ) = fl(x q ). 

Clearly, p > and //(x^ 1 ) / //(a*) = s. So if s = 1 then /^a*" 1 ) < x* < /i(a*) 
and we deduce that G(F) has a positive arc from j to i. Similarly, if s = —1 then 
fi{x p ~ v ) > Xj > fi(x p ) and we deduce that has a negative arc from j to i 
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Case j = i. By definition, s = f' i {x)f' i (Fi(x)) and f[{x) / f'^F^x)) thus s = — 1. Suppose 
that f'i{x) > 0, the other case being similar. Then q = fi(x) — X{ > and //(Fj(x)) < 0. 
Since x° = x and x q = Fi(x), we deduce that x® < fi(x°) = x\ and fi(x q ) < x\. Thus, 
there exists a smallest < p < q such that fi{x p ) < x\. Clearly, p > and x\ < fi(x p ^ 1 ). 
Thus fi{x p ) < fi{x p ~ l ) and we deduce that G{F) has a negative arc from i to itself. □ 

Lemma 2 Let x 1 , . . . , x r ) be an elementary path of T(F) of length r > 1, and let 
i £ If{x t ). If f' i (x p ) 7^ f'i{x r ) for all < p < r, then there exists j € If{x°) such that 
|JgZo^F(^ 9 ) ho-s a path from j to i of sign /j(x°)//(x r ). 

Proof — We reason by induction on the length r of the path. 

Case r = 1. Since x 1 ) is an arc of T(F) there exists j € If{x°) such that x 1 = Fj(x°). 
Following the conditions of the lemma //(x°) / fifa 1 ), and thus, by definition, ^(x°) has 
an arc from j to i of sign /j(x°)/ J '(x 1 ). 

Case r > 1. Since (x r_1 ,x r ) is a path of of length 1 satisfying the conditions of 

the lemma for i € If(x v ), following the base case, there exists k G Ii?(x r_1 ) such that 
^F(x r_1 ) has a path from k to i of sign 

Now, consider the smallest < p < r such that f' k {x p ) = /(,(x r-1 )- First, suppose that 
p = 0. Then k € If{x°) and f' k ( K x )fl(x r ) is equals to sign Sfcj of the path of ^i?(x r ^ 1 ) from 
k to i mentioned above, so that the lemma holds. Now, suppose that p > 0. Then, by the 
choice of p, for all < I < p, we have f' k {x l ) ^ f'k{x p )- Thus, the path (x°, . . . ,x p ) satisfies 
the conditions of the lemma for k G If{x p ). Since p < r, by induction hypothesis, there 
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exists j £ If{%°) such that U«=o^H x9 ) has a path from j to k of sign 



Since ( x 



.r-l 



) contains a path from to i of sign s^j, we deduce that Ug=o ^f(^ <7 ) contains 



a path from j to i of sign 



and since f' k (x p ) = f' k (x r 1 ), we deduce that Sji = f'Ax°)f' i (x r ). 



□ 



Lemma 3 Let A be a cyclic attractor ofT{F). If there exists x G A such that \Ip(x)\ = 1 



Proof — Suppose that there exists x° € A such that = 1, and let i be the unique 

element of Ip(x ). Suppose that fl(x°) > 0, the other case being similar. Let x 1 = Fi(x). 
Then T(F) has an arc from x° to x 1 and we have x® < x\. Since x° S A, we have x 1 £ A, 
and we deduce that T(F) has an elementary path ( X , X • • • • • X r ) from x 1 to x r = £° whose 
all the vertices belong to A. If fl(x p ) > for all < p < r, then < x p+1 for all < p < r, 
and we deduce that xj < x\ = x®, a contradiction. Thus, there exists a smallest < p < r 
such that fl(x p ) < 0. Then, (ar, ar, . . . , ar 1 ) is an elementary path where i E Ip(x p ) and by 
the choice of p, we have //(ar) 7^ fl(x p ) for all < Z < p. So, according to Lemma [21 there 
exists j € If(x°) such that Ug=o^H x<? ) contains a path from j to i of sign fj{x°)f' i (x p ). 
Since If^ ) = {«}> we have j = i and consequently, \JiqZo^F(x q ) contains a path from i 
to itself, and thus a circuit, of sign f' i {x°)f' i {x p ). By construction, fl(x°)fl(x p ) < 0, thus 
this circuit is negative, and since . . . , ar 1 " 1 } C A, it is contained in (J xgj4 ?fp(x). □ 
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£/ien Ukga &f(x) has a negative circuit. 



Lemma 4 Let A be a cyclic attractor ofT(F). If \Ip(x)\ > 1 for all x £ A, then there 
exists H : X — > X such that T(H) contains a cyclic attractor strictly included in A, and 
such that &h(x) is a subgraph of^p(x) for all x £ X. 

Proof — Suppose A to be a cyclic attractor of T(F) such that \Ip(x)\ > 1 for all x £ A. 
Let y be any state of A. Then If{v) contains at least two elements, and without loss of 
generality, we can suppose that 1 £ If{v)- Consider the map H : X — > X defined by: 

Vx G X, H(x) = (h 1 (x),h 2 (x),...,h n (x)) = (xi, / 2 (x), . . . , /„(x)). 

We first prove that A is a trap domain of T(H). For that, it is sufficient to prove that, 
given any x £ A and i G Ih(x), we have Hi(x) G A. Since h\(x) = x±, 1 Ih(x), so i 7^ 1. 
Thus Fj(ar) = Hi(x), and since ^4 is a trap domain of we have Fj(ar) G A and we 

deduce that Hi(x) G i as expected. So A is a trap domain of T(H) and, by definition, 
T(H) contains at least one attractor B C A. 

We claim that B is a cyclic attractor of T(H). Let x £ B. Then x G A so > 1 

and we deduce that Ip(x) contains an index i 7^ 1. Then, Xj 7^ /i(x) = /i«(x) sox/ Hi(x). 
Since x £ B we have Hi(x) G -B. So |J3| > 2, i.e. 5 is a cyclic attractor of T(H). 

We now prove that B C A (strict inclusion). Suppose, by contradiction, that B = A. 
Since 1 € /^(y) and !/ £ 4, we have y ^ F\{y) £ A = B. Since I? is an attractor of 
r(iJ), we deduce that r(il) has a path (x^x 1 , . . . ,x r ) from x° = y to x r = F\(y). Since 
/ii(x) = xi for all x G X, we have x? = x\ = ■ • • = x\. So yi = fi(y), a contradiction. 

It remains to prove that £f#(x) is a subgraph of ^p(x) for all x G X. If (j,s,i) is an 
arc of Sf#(x), then by definition, foj(x) 7^ and h' i {Hj(x)) 7^ 0. So j 7^ 1 and i 7^ 1. Thus 
/j = hj and /j = /ij. It is then clear that (i, s,j) is an arc of Gf(x). □ 
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Lemma 5 If A is a cyclic attractor ofT(F), then UxgA^p( x ) has a negative circuit. 

Proof — Let U be the set of couples (F, A) such that F is a map from X to itself, and such 
that A is a cyclic attractor of T(F). Let -< be the well funded strict order on U defined by 
(H, B) ~i (F, A) if and only if B is strictly included in A. Reasoning by induction on the 
set U ordered by -<, we show that, for all (F, A) 6 U, \J x£A ^f(x) has a negative circuit. 

Base case. Let (F, A) be a minimal element of (U, -<). If > 1 for all x E A, then, 

following Lemma HI there exists (H, B) G U such that (H, B) -< (F, A), and this contradict 
the minimality of (F, A). So there exists x G A such that = 1 and, following 

Lemma [3 LbeA^H^) nas a negative circuit. 

Induction step. Let (F, A) be a non-minimal element of (U, -<). By induction hypothesis, 
for all (H,B) -< (F,A), \J x€B &h(x) has a negative circuit. If, for all x € A, we have 
> 1, then following Lemma HJ there exists (H,B) -< (F, yl) such that &h( x ) is a 
subgraph of ^f(^) for all ifl. Since B C A, we deduce that [j xeB ^h{x) is a subgraph 
of U a ;GA^ 7, ( x )' ana - since, by induction hypothesis, \J x€ b^h(x) nas a negative circuit, we 
deduce that \J x£A &f(x) has a negative circuit. Otherwise, there exists x & A such that 
| Ji^(a:) | = 1, and following Lemma[3l UxeA^H 3 ') has again a negative circuit. □ 

Proof of Theorem 1 — If A is a cyclic attractor of then by Lemma[5j \J x& a^f(x) 

has a negative circuit. By Lemma [IJ UzgA^H x ) is a subgraph of G(F) and we deduce 
that G(F) has a negative circuit. □ 

Remark 7 The key lemma is clearly Lemma 5, which shows that it is sufficient to consider 
the restriction of F to a cyclic attractor A in order to obtain a negative circuit. 
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4 A variant for gene regulatory networks 

In this section, we establish a variant of Theorem 1 that is more suited to the modeling 
of gene networks. To model the behaviors of a network of n genes, Thomas [26] [28] 
proposes to consider an "unitary" asynchronous state transition graph T[F] that is slightly 
different than T(F). In T[F], each transition starting from a given state x involves, as 
in r(F), the evolution of the state Xi of at most one component i € If(x), but in F[F], 
this state Xi is not updated to fi(x): it is increased or decreased by a unit depending on 
whether X{ < fi{x) or x« > fi(x). Thanks to this updating rule, unitary asynchronous state 
transition graphs can be seen as discretizations of piece-wise linear differential equation 
systems [151 120]. 

Definition 6 The unitary asynchronous state transition graph of F, denoted T[F], is the 
asynchronous state transition graph T(F) of the map F : X — > X defined by 

F(x) = (fi(x),...,f n (x)), fi(x) = Xi + f-(x) (i = l,...,n). 

Remark 8 In the Boolean case, T[F] = T(F). 

We are now confronted to the following problem: G(F) cannot be seen as the interaction 
graph of the network whose dynamics is described by T[F], since maps H such that G(H) ^ 
G(F) and T[H] = T[F] may exist. In addition, it is not satisfactory to see G(F) as the 
interaction graph of the network whose dynamics is described by T[F], since maps H such 
that G(H) is a strict subgraph G(F) and such that T[H] = T[F] may also exist. 

To solve this problem, Richard and Comet [T3] define a subgraph G[F] of G(F) that only 
depends on T[F] and provide, in this way, a natural and non-ambiguous definition of the 
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interaction graph of the network whose dynamics is described by T[F]. Furthermore, one 
can show that G[F] is, with respect to the subgraph relation, the smallest interaction graph 
from which one can obtain T[F] by following the logical method developed by Thomas to 
model gene networks |12| , 

Definition 7 We denote by G[F] the interaction graph that contains a positive arc from 
j to i if there exists x G X with Xj + 1 G Xj such that 



and that contains a negative arc from j to i if there exists x G X with Xj + 1 G Xj such that 



Remark 9 G[F] is a subgraph of G(F), and in the Boolean case, G[F] = G(F). 

We now establish, in this setting, the following discrete version of second Thomas' 
conjecture (which is, as Theorem 1, an immediate consequence of Lemma|5|: 

Theorem 2 IfT[F] has a cyclic attractor, then G[F] has a negative circuit. 

Lemma 6 For all x G X, &p(x) is a subgraph of G[F\. 

Proof — First observe that f[(x) = f[(x) for all x G X and i G {1, . . . , n}. Furthermore, 
if fi{x) < Xi (resp. fi(x) > xt) then fi(x) < fi(x) (resp. fi(x) > fi(x)). 

Now, suppose that &p{x) has an arc from j to i of sign s with j ^ i. Let 



fi(xi,-.. 



, Xj, 



X n ) ^ Xi < fi (xi , . . . , Xj + 1, . . . , Xn) , 



fi{x x ,... 



, Xj , • • • , X n ) > Xi ^ fi (X\ , . . . j Xj + 1 , . . . , x n 



y = Oi,. . . 




) 
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and observe that y = Fj(x). Suppose that f[{y) > 0, the other case being similar. Then, 
by definition, fj(x) = s and f[{x) < 0. Thus fi(x) < x« = yi < fi(y) and we deduce that 

fi(x) < fi( x ) <xi = yi< My) < fi{y)- 

So if fj(x) = s is positive then 

fi(x) <Xi< fi(y) = fi(x!,.. .,Xj + l,...,x n ) 
and we deduce that G[F] has a positive arc from j to i, and if fj(x) = s is negative then 

fi(yi, • • • , Vj + l, • • • , 2/n) = /i(a;) < ^ < /i(y) 

and we deduce that G[F] has a negative edge from j to i. 

Suppose now that 2fp(x) has an arc from i to itself of sign s. By definition, we have 
s = f' i {x)f' i (Fi(x)) and fl(x) / fl(Fi(x)) so that s is negative. Suppose that fl(x) > 0, the 
other case being similar. Then, Fi(x) = (xi, . . . , X\ + 1, . . . , x n ) and //(Fj(x)) < 0. Thus 

fii.Xli ■ ■ ■ -i Xi + 1, . . . , X n ) < Xj < /j(x) 

and we deduce that 

/i(xi,. . . ,Xi + 1,. . . ,X n ) < /i(xi, . . . ,Xi + 1, . . . ,X n ) < Xi < /j(x) < /j(x). 

Consequently, has a negative arc from i to itself. □ 
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Proof of Theorem 2 - Since T[F] = T(F), if T[F] has a cyclic attr actor A, then by 
Lemma 5, Uxex^F'( a ') nas a negative circuit. Following the previous lemma, LLex^M^) 
is a subgraph of G[F], and we deduce that G[F] has a negative circuit. □ 

Corollary 2 If G[F] has a no negative circuit, then F has at least one fixed point. 

Proof - If F has no fixed point, then T[F] has at least one cyclic attractor, and following 
Theorem 2, G[F] has a negative circuit. □ 

Remark 10 Since G[F] is a subgraph of G(F), Corollary 2 is stronger than Corollary 1 
(the same conclusion is obtained under a weaker condition). In addition, from Theorems 1 
and 2, it is clear that: if T(F) or T[F] has a cyclic attractor, then G(F) has a negative 
circuit. This generalizes Theorem 1 (the same conclusion is obtained under a weaker 
condition). Indeed, as showed by the following two examples, the presence of a cyclic 
attractor in T(F) (T[F]) does not imply the presence of a cyclic attractor in T[F] (T(F)). 

Example 2 n = I, X = {0, 1, 2} and F defined by F(0) = 2, F(l) = 1 and F{2) = 0. 
The state transitions graphs T(F) and T(F) are the following: 

F(F) T[F] 



> 1 «e 



We see that T(F) has a cyclic attractor and that T[F] has no cyclic attractor. The interac- 
tion graph G(F) is the interaction graph with one vertex and a negative arc from this vertex 
to itself: it has thus a negative circuit. The interaction graph G[F] is the interaction graph 
with one vertex and no arc (it is a strict subgraph ofG(F)). This shows that the presence 
of a cyclic attractor in T(F) does not imply the presence of a negative circuit in G[F]. 
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Example 3n=l,I = {0,l, 2} and F defined by F(0) = 0, F(l) = 2 and F(2) = 0. 
The state transitions graphs T(F) and T(F) are the following: 




We see that T[F] has a cyclic attractor and that T(F) has no cyclic attractor. The inter- 
action graphs G(F) and G[F] are equal to the interaction graph with one vertex and both a 
positive and a negative arc from this vertex to itself (G{F) and G[F] have thus a negative 
circuit). 

5 Concluding remarks 

The weakest condition allowing the asynchronous iterations of F to describe sustained 
oscillations is the presence of a directed cycle in T(F). However, as showed by the following 
example, the presence of a directed cycle in T(F) does not imply the presence of a negative 
circuit in G(F) (one can only show that it implies the presence of a circuit in G(F)). This 
shows that structures in T(F) stronger than directed cycles (such as cyclic attractors) are 
needed to obtain a negative circuit. 
Example 4 n = 3, X = {0, l} 3 and F is defined by 

fi{x) = x 3 
f 2 {x) = Xl 

h{x) = x 2 . 



19 



The asynchronous state transition graph T(F) (which is here equal to T[F]) and the inter- 
action graph G(F) (which is here equal to G[F]) are the following: 



F(F) 



(0,1, 0)« 



(1,1,0) 



\ 




G(F) 



(0,1,1) 



>(1,1,1) 



i 



® 



(0,0,1) 



>(i,o,i) 




\ 



(0,0,0)< 



(1,0,0) 



We see that T(F) has a directed cycle and that G(F) has no negative circuit. 

A second remark is that it is not easy to find other classes of iterations for which 
Theorem 1 remains valid. Consider for instance the synchronous state transition graph 
A(F) that encodes the behaviors of the iteration x t+l = F(x t ): the set of vertices of A(F) 
is X and the set of its arcs is {(x, F(x)) \ x 6 X, x ^ F(x)}. The cyclic attractors of such a 
(deterministic) state transition graph A(F) are naturally defined to be the directed cycles 
of A(i< 1 ). However, the following example shows that the presence of a directed cycle in 
A(F) does not imply the presence of a negative circuit in G(F) (Robert [T7] prove that 
it only implies the presence of a circuit in G(F)). 
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Example 5 n = 2, X = {0, l} 2 and F is defined by 



fi{x) = x 2 
f 2 {x) = Xl. 

The synchronous state transition graph A(F) and the interaction graph G(F) are as follows: 

A(F) G(F) 




We see that A(F) has a cyclic attractor and that G(F) has no negative circuit. 

Finally, we can ask if, under the condition that T(F) has a cyclic attractor, a conclusion 
stronger than U G(F) has a negative circuit" could be obtained. Following Example 2, the 
presence of a cyclic attractor in T(F) does not imply the presence of a negative circuit 
in the subgraph G[F] of G(F). So, another direction has to be taken. As showed below, 
previous results on the links between the interaction graph and the dynamical properties 
of automata networks suggest to improve the conclusion of Theorem 1 by studying if the 
presence of a cyclic attractor in T(F) implies the presence of a negative circuit in a local 
interaction graph associated with F. 

Definition 8 For all x <G X, the local interaction graph of F evaluated at state x is the 
interaction graph Gf{x) that contains a positive (negative) arc from j to i ifxj + 1 G Xj and 

fi(x±, . . . , Xj + 1, . . . , X n ) fi{x\, . . . , Xj, . . . , X n ) 
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is positive (negative), or if Xj — 1 G Xj and 

fi \Xl j • • • > •Ej j • • • 3 X n ) fi[.Xl ) ■ ■ • j •Ej 1 3 • • • > -^n ) 

is positive (negative). 

Remark 11 Gf(x) is a subgraph ofG(F). More precisely, G(F) = (J xg jf Gf(x). 

With this material, Richard and Comet [13] prove the following local version of first 
Thomas' conjecture: 

Theorem 3 [13J IfT[F] has several attractors, and in particular if F has several fixed 
points, then there exists x € X such that Gf(x) has a positive circuit. 

Let us also mention the following fixed point theorem proved by Richard [13] (and previ- 
ously proved by Shih and Dong [18] in the Boolean case): 

Theorem 4 |14j If Gp{x) has no circuit for all x € X , then F has a unique fixed point. 

The proof of Theorem 4 done in [J5] reveals that if Gp(x) has no circuit for all x E X, 
then F has a unique fixed point £, and, in addition, for all T[F] has a path from x 

to £. It is then clear that the presence of a cyclic attractor in T[F] implies the presence of 
a circuit in Gf{x) for at least one x € X. We then arrive to the following natural question: 

Question 1 Does the presence of a cyclic attractor in T[F] or T(F) implies the presence 
of a negative circuit in Gf{x) for at least one x € X? 

Clearly, a positive answer would improve significantly Theorem 1 or 2 by providing a local 
version of second Thomas' conjecture. However, the following example shows that the 
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answer is negative. This highlights the fact that it is necessary to take a union of local 
interaction graphs in order to obtain, from a cyclic attractor, a negative circuit. 

Example 6 n = 2, X = {0, 1, 2, 3} 2 and F is defined by: 



3 if X2 = 3 or if X2 > and x\>2 
otherwise 



h{x) = { 



3 if Xi = or if x\ < 3 and X2 > 2 
otherwise 



The asynchronous state transition graph T(F) is the following: 




(2,3) >(3,3) 

A 



(2,2) >(3,2) 



(2,1) >(3,1 



Y V V 

(0,0) < (1,0) (2,0) (3,0) 
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The unitary asynchronous state transition graph F [F] is the following: 



(0,3) > (1,3) > (2,3) > (3,3) 

AAA 

\|/ 

(0,2) < (1,2) (2,2) >(3,2) 

A 

(0,1) < (1,1) (2,1) >(3,1) 

A 

\|/ \|/ \|/ 
(0,0) < (1,0) < (2,0) < (3,0) 



The interaction graph G{F), which is here equal to G\F\, is the following: 



:© 



We see that {(0, 0), (0, 3), (3, 3), (3, 0)} is a cyclic attractor ofF(F) and that G(F) has a neg- 
ative circuit. We see also that {(0, 0), (0, 1), (0, 2), (0, 3), (1, 3), (2, 3), (3, 3), (3, 2), (3, 1)(3, 0)} 
is a cyclic attractor ofF[F] and that G[F] has a negative circuit. However, for all x € X, 
the local interaction graph Gf(x) has no negative circuit. Indeed, forx G {(1, 0), (0, 0), (0, 1)} 
and x G {(2, 3), (3, 3), (3, 2)} ; Gf{x) is as follows: 
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for x £ {(3, 1), (3, 0), (2, 0)} and x £ {(0, 2), (0, 3), (1, 3)}, G F (x) is as follows : 



i 



/or a; = (1, 1) and x = (2, 2), Gf{x) is as follows: 



O 



2 i 



/or x = (1, 2) and x = (2, 1), Gf (x) zs as follows: 



*0 



The fact that Theorem 3 establishes the uniqueness of a fixed point for F under the condi- 
tion that Gf(x) has no positive circuit suggests the following weaker version of Question 1: 

Question 2 Does the absence of a negative circuit in Gf{x) for all x € X implies the 
presence of at least one fixed point for F ? 

A positive answer would improve significantly Corollary 1, and would give, together with 
Theorem 3, a very nice "dichotomous" proof of Theorem 4. However, the previous example 
shows that Question 2 as also a negative answer. Nevertheless, 



Questions 1 and 2 remain open in the Boolean case. 
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